Spin-polarized stable phases of the 2-D electron fluid at finite temperatures. 
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The Helmholtz free energy F of the interacting 2-D electron fluid is calculated nonperturbatively 
using a mapping of the quantum fluid to a classical Coulomb fluid [Phys. Rev. Letters, 87, 206 
(2002)]. For density parameters r 3 such that ~ 25 < r s , the fluid is unpolarized at all temperatures 
t = T/Ef, where Ef is the Fermi energy. For lower densities, the system becomes fully spin 
polarized for t smaller than ~ 0.35, and partially polarized for ~ 0.35 < t < 2, depending on the 
density. At r s ~ 25 — 30, and t ~ 0.35, an "ambispin" phase where F is almost independent of 
the spin polarization is found. These results support recent claims, based on quantum Monte Carlo 
results, for a stable, fully spin-polarized phase at T = for r s larger than ~ 25 — 26. 
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Introduction. The physics of the uniform two- 
dimensional electron gas (2DEG) depends crucially on 
the "coupling parameter" T = (potential energy) /(kinetic 
energy) arising from the Coulomb interactions. The T 
for the 2DEG at T = and mean density n is equal 
to the mean-disk radius r s = (im)^ 1 ^ 2 per electron. In 
the absence of a magnetic field or disorder effects, the 
density parameter r s , the spin polarization £ and the 
temperature T are the only relevant physical variables. 
As the coupling constant increases, the 2DEG becomes 
more like a "liquid" , but here we continue to loosely call 
it an electron "gas". A knowledge of the phase dia- 
gram at arbitrary r s , £ and T is of crucial importance 
to a clear understanding of many intriguing phenomena 
(metal-insulator transition^, 2-D electrons in CuO lay- 
ers, spin-dependent processes etc.) associated with the 
2DEG. 

A number of recent studies^, ||, ||| have suggested 
that the 2DEG at T =0 becomes fully spin-polarized 
for r s > r c with r c ~ 25 — 30. However the energy 
difference between the unpolarized and polarized phases 
is very close to the statistical uncertainty of the meth- 
ods used. We have completed a study of this transition 
at finite temperatures where the energy differences for 
a larger data base become more significant and reveal 
a rich phase diagram which is fully consistent with the 
claimed T — polarized-fluid regime. It turns out that 
the partially degenerate strongly correlated 2DEG sup- 
ports partially spin-polarized phases. This is particularly 
interesting since the 3DEG has also been recently shown 
to have partially spin-polarized stable phases^] 

The finite- T phase diagram is examined using a com- 
putationally simple, conceptually novel method for cal- 
culating the quantum-PDFs and other properties (e.g, 
static response). The method has been applied to the 
2DEG§, the 3DEG[§ 0, and hydrogen fluids @, at ar- 
bitrary coupling, spin polarization and temperature. It 
is based on identifying a "quantum temperature" T q such 
that the correlation energy of the corresponding classi- 
cal Coulomb fluid at T q is equal to that of the quantum 



fluid at T = 0. After that, hypernetted-chain (HNC) 
methods are used for the classical fluid. This classical 
mapping of quantum fluids within the HNC was named 
the CHNC. The application of the CHNC method to the 
2DEG required the inclusion of short-range correlation 
effects ("bridge terms") going beyond the usual HNC 
approximation. The resulting CHNC procedure for the 
2DEG proved to be remarkably accurate In fact even 
an application of CHNC at T = 0, and without bridge 
corrections was examined by Bulutay and Tanatar and 
found to be very useful 0. 

We use Hartree atomic units, and consider a unit area 
A with n particles per unit area, with n = l/(irr 2 ). The 
number density of each spin species is n.; = n(l + cr^)/2, 
where 1=1,2 would indicate up or, down spins with 0{ = 
±1. The spin polarization £ = (rii — n^jn. The total 
Helmholtz free energy per unit area, F, has the form 
Fa + F x + F c , where Fq, F x , F c are the non-interacting, 
first-order exchange, and correlation, contributions to F. 
Thus all effects beyond first-order exchange are lumped 
into F c , a convention used in density functional theory 
(DFT). The DFT exchange-correlation free energy would 
hence be F xc — F x + F c . In indicating spin-components 
we use the notation F°,F^, and Ff etc. The canonical 
ensemble is used so that the Helmholtz free energy is 
directly determined, unlike in the approach where the 
grand-potential Q is determined in terms of the chemical 
potential which has to be inverted to obtain the density. 
This means Fq and F x (being the first-order exchange) 
are determined in terms of the non-interacting chemical 
potential fj,o, while F c has to be determined via a coupling 
constant integration, as discussed by us in r of § [Uj| and 
in Dundrea et alfll2[. 

The non-interacting free energy. Isihara and others 
have considered the expansion of the free energy about 
T = and obtained the leading terms in r s , and T 
]l0| . However, the thermodynamic functions near T = 
have a particularly unfavourable structure for construct- 
ing such expansions. Since our calculations are nonper- 
turbative, we evaluate Fq etc., numerically, without ex- 
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pansion. 

Using the the relation, F = E — TSq, where E and 
So are the noninteracting internal energy and the en- 
tropy, Fq can be calculated. An independent evaluation 
is obtained from the relation Fq — njj,a ~ Eq. Also, at 
T = 0, we have, in Hartrees 



F°/n = E°/n=(l + a. l C) 2 /4r 2 s 



(1) 



Defining the Fermi function /(e) = 1/(1 + exp(J3e — rj)) 
where (3= 1/T, rj = fjP/T, 



m = (T/27r)7 (r ?i ) 
E t = (T 2 /pi)h( Vl ) 

where the Fermi integral I v (z) is defined as 

. dxx" 

h{z) 



(2) 
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la 1 + e x - z 

The Iq integral can be evaluated explicitly to give: 

Vl = log[exp(l/ti) - 1], U = Ef/T (6) 
Ef = E F (l + ai(), E F =irn (7) 

(8) 

The 1\ integral can be evaluated directly or reexpressed 
in terms of the dilogarithm function. We have evaluated 
Fa by several methods; the dilog function was evaluated 
using the CERN library routine. 

The exchange contribution to the free energy. The first- 
order unscreened exchange free energy F x consists of F x , 
where i denotes the spin species. At T = these reduce 
to the exchange energies: 



Ef/n = 

Hence, 

E x /n=(Ef +E%)/n 
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where x\ and X2 are the fractional compositions (l±£)/2 
of the two spin species. At finite T, we have: 



3 q/9 /"" I 2 U9 (u)du 
F*/E? = —t 3 / 2 1 1 - 
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(11) 



Here /_i/2 is the Fermi integral. The total exchange free 
energy is F x = SFf. The accurate numerical evaluation 
of Eq. [ll] requires the removal of the square-root singular- 
ity by adding and subtracting, e.g, I 2 (— \r)\)/(v— \eta\y/ 2 
for the case where i] is negative, and v = u, and so on. 

A real-space formulation of F x = Ff + using the 
zeroth-order PDFs fits naturally with the CHNC ap- 
proach. Thus 



F x /n = n I ^YhUr) 



r / j v ' 

i<j 



(12) 



Here /i°-(r) = g^ir) — 1. In the non- interacting system 
at temperature T, the antiparallel h\ 2 i y i z -i 3i2( r i T) — 1, 
is zero while 



1 



h° n (v) = ,E kl , k! n(fc)n(fe)e ,(kl 



■lf(r)f 



Here k,r are 2-D vectors and n(k) is the Fermi oc- 
cupation number at the temperature T. At T = 
f(r) = 2Ji(kir)/k r where J\{x) is a Bessel function. We 
have numerically evaluated the exchange free energy by 
both methods, i.e, Eq. [ll] and Eq. [l2|, as a numerical check. 

Correlation free energy. The correlation contribution, 
F c is evaluated together with the exchange contribution 
F x via a coupling constant integration over the distribu- 
tion functions as follows: 



F xc /n = J dXn J 



2nrdr 



h^X) (13) 



Thus we need the interacting PDFs gij (r, A) for each r s , 
C and T for a number of values (we have used 7 values) of 
the coupling constant A. The PDFs were evaluated using 
the CHNC method, as described in Ref.[Q. We present 
a brief summary of the CHNC method in the context of 
the present study. 

Brief summary of the CHNC method. The essence of 
the CHNC method is to start from the quantum mechan- 
ical g (r) of the non-interacting problem and build up the 
interacting g(r) by classical methods. 

The fluid of mean density n contains two spin species 
with concentrations X{ — rii/n. While T is the physical 
temperature of the 2DEG, we consider a classical fluid of 
temperature T c f = 1//3 C /- Since the leading dependence 
of the energy on temperature is quadratic, we assume 
that T c f= a/(T 2 + Tjj) where T q is the quantum correc- 
tion. This is clearly valid for T = and for high T, 
and was justified in more detail in ref. The "quantum 
temperature" T q was shown to be given 



t = Tq/Tp = 2/[l + 0.86413(r 
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where the Fermi temperature Tp is simply the Fermi en- 
ergy Ep = 1/r 2 in our units. In effect T q is really a 
single-parameter representation of the density-functional 
2D correlation energy. 

The HNC equations for the PDFs of a classical fluid, 
and the Ornstein-Zernike(OZ) relations are|lg|]: 

gij(r) = exp[-p c f<f>ij(r) + hij(r)-Cij(r)+Bij(r)] 
hij{r) = c i3 (r) + Y, s n s / dr'h ilS (\r - r'\)c s j{r') (15) 



This equation involves two quantities, viz., (i) the pair- 
potential 4>ij{r), (ii) the bridge function i?^- (r) |l4|, p^[ . 
The other terms, e.g, c(r), is the "direct correlation func- 
tion" of the OZ equations, while hij(r) = <7y (r)-l. 

Pair-potential. The <f>ij (r) is the pair potential between 
the species i, j . For two electrons this is just the Coulomb 
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potential V cou (r). If the spins are parallel, the Pauli 
principle prevents occupation of the same spatial orbital. 
This effect is equivalent to a "Pauli exclusion potential" , 
V(r) which when used in the HNC equation reproduces 
the g%(r). Thus <£ y (r) becomes V{r)8 io + V cou (r). The 
Coulomb potential V CO u (r) for a pair of point-charge elec- 
trons is 1/r. However, an electron at the temperature T c f 
is localized to within a thermal wavelength. Thus, we use 
a "diffraction corrected" form 



V cou (r) = (l/r)[l-e 



-rk th 



(16) 



In the case of the 2DEG kth was explicitcly determined by 
numerically solving the Schrodinger equation for a pair 
of 2-D electrons in the potential 1/r and calculating the 
electron density in each normalized state That is, by 
solving the 2-D Schrodinger equation (only s -waves are 
needed for the r — > case): 

d 2 R IdR n fc 2 1. n 

l? + rdF + 2m (T-J R = (17) 

Here m* = 1/2 is the effective mass of the electron pair. 
The radial function has the asymptotic form: 

R(r, k) - [2-]V*coa(kr _ H£l„( r ) + 8 2d ) (18) 
irkr k 

The "on-top" density n(0) for electrons at the tempera- 
ture 1/(3 is: 



n(0) _ J °° k 2 dkexp{-/3k 2 /2m*)R 2 (0,k) 
~n~ ~ J °° k d dkexp{-(3k 2 /2m*) 



(19) 



Using the same form for the diffraction correction, i.e., 
Vrf(r)=exp(— kthf)/fj h was found that 



hh/k° tl 



th 



1.158T C / 03 



where T c f is in a.u., and k® h is the de Broglie thermal 
momentum (27rm*T c ^) 1 / 2 . 

Bridge contribution. The "bridge" term arises from 
higher-order cluster interactions. They seem to play 
a role similar to the "back-flow corrections" used in 
QMC trial wavefunctions Jl7| . Since the Pauli exclusion 
already acts to reduce clustering of parallel spin elec- 
trons, the Ba bridge corrections are unimportant. In 
Ref. Q we included the Bn terms via a simple hard- 
disk bridge function, jl5j with a hard-disk packing frac- 
tion ri p = 0.1175r s 1/3 (T c/ /£ F ). Since t = T cf /E F is itself 
a function of r s , the zero temperature packing fraction 
i] p (T — 0) is just a function of r s . In generalizing the 
above expression to finite T, we note that r s is simply 
the coupling constant 7 = P.E/K.E per electron. We 
take the potential energy P.E = l/r s for a pair of elec- 
trons, and evaluate the kinetic energy (K.E.) from the 
non-interacting internal energy as in Eq. §. Then the fi- 
nite temperature hard-disk packing fraction is given by 




T) p = 0.235 7 1/3 /[l + 0.86413( 7 1/e - 1) 



1/6. 



(20) 



FIG. 1: The phase diagram of the 2DEG as a function of 
the density parameter r s and temperature, showing fully and 
partially spin-polarized stable phases. The region near the 
dotted lines near r 3 ~ 25 and t ~ 0.35 represents an "am- 
bispin" phase where the free energy is almost independent of 
the spin polarization. 



The next step is to use the faj (r) , and solve the coupled 
HNC+Bridge and OZ equations for the two-component 
interacting fluid. Evaluation at several values of the cou- 
pling constant and integration as in Eq. [l3] leads to the 
correlation free energy. 

Results. The fluid-phase with the lowest Helmholtz 
free energy F(r s ,T, () is the stable phase. We have eval- 
uated F for the range t =T/E F < 2, and 5 < r s < 40, 
and usually for five values of £ and found that the sta- 
ble phase is the unpolarized (paramagnetic) phase for 
higher temperatures and densities. The initial sugges- 
tion of Varsano and Senatore that the 2DEG at T=0 is 
spin polarized near r s ~ 30, and the more recent dis- 
cussion by Attaccalite et al., that a transition occurs at 
r s =25.7 are greatly strengthened by our results at finite- 
T. At T=0 we find a ( = ^ ( = 1 transition at 
r s = r c (( : — ► 1) = 26.7 ± 0.8, but the stabilization 
energy is within the inherent uncertainities of the CHNC 
method. However, when taken together with all the re- 
sults for neighbouring finite-T values, a more plausible 
picture emerges ( Fig. [l]). When the temperature is in- 
creased, the value of r c (£ : — > 1) decreases somewhat 
and then begins to increase. However, when t reaches 
~ 0.35, partially polarized states with £ ~ 0.7 to 0.9 
begin to compete for stability. This partially polarized 
state remains stable as long as the temperature is suffi- 
ciently low. An interesting aspect of the partially polar- 
ized regime is the existence of "ambi-spin" regions where 
the spin projection is NOT unambiguously determined 
since the energy differences as a function of the polar- 
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FIG. 2: The polarization dependence of the Helmholtz free 
energy plotted as AT((") = F(r 3 ,t,() — F(r 3 ,t,0) at several 
temperatures and densities close to a tranition (see the text). 



ization are very small ( see below). For t >^ 1.5, the 
stable phase becomes the usual unpolarized (paramag- 
netic) phase. 

A plot of the energy difference AF(r s ,t, () — T(£) — 
T(0), as a function of the spin polarization is shown in 
Fig. |^ at several temperatures and r s values close to one 
of the transition lines. Thus the t — 0.05 curve is at r s ~ 
25.86 which corresponds to r c (£ : — > 1). This curve is 
similar to that given in Fig. 1 of Ref. p). However, the 
the transition regions near r s ~ 26 and t ~ 0.35 provides 
an example of a phase where the Helmholtz fee energy 
is almost independent of the spin polarization. Such an 
"ambispin phase" would be very sensitive to any external 
effects (impurities etc.,) that we have not considered in 



this study. The regions in Fig. |l] close to the dotted 
line, begining from the "triangular" region near r s ~ 
25 and t ~ 0.35 are such regions. The partially spin 
polarized state is quite unusual since, in some regions, a 
very light deviation from £ = 1 leads to a significantly 
large stability. Thus the curve at t = 0.5 is typical of 
densities near r s ~ 30 — 32 where the partially polarized 
stable is very close to unity but differs slightly from unity. 

At T = the liquid phase becomes a Wigner crystal 
near r s ~ 35. However, as the temperature is increased, 
this value will also increase. At low temperatures the 
crystal would be in equilibrium with the fully spin po- 
larized phase. We have not investigated the equilibrium 
between the solid and liquid phases. 

Thus the finite-T 2DEG system provides many intrigu- 
ing questions while complementing our understanding of 
the T = system. Our computer codes for generating 
the PDFs, the f xc ((r s ,T, () and related properties of the 
2D electron fluid may prove to be useful for further stud- 
ies of this system. The codes may be accessed via the 
internet @. 

In conclusion, we have evaluated the total Helmholtz 
free energy F(r s ,T, () — F + F x + F c , in the canonical 
ensemble, for a range of r s < 40, t < 2, and < z < 1. 
The terms Fq and F x were evaluated by several methods, 
without resorting to perturbation expansions, while the 
F c was evaluated via a coupling constant integration of 
the pair-distribution functions obtained from the CHNC 
equations. The finite-T results are consistent with recent 
claims, based on quantum Monte Carlo calculations, of 
the existence of a fully spin-polarized stable phase at T 
=0 and r s ^ 25 — 26. Our results show that the finite-T 
phase diagram of the clean, uniform 2DEG contains sta- 
ble fully and partially spin-polarized phases. Ambispin 
phases where the free energy is quasi-independent of the 
spin polarization arc also found. It is hoped that these re- 
sults would stimulate new experiments and finite-T QMC 
simulations. 
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